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Imaging 

This invention relates to imaging. 

5 

The invention is concerned especially with imagining 
techniques using CT (computed tomography) , otherwise 
known as CAT (computer assisted tomography) , scanning. 

10 CT scanning provides an image that may be displayed on a 
monitor, of a transverse slice of the head or body of a 
patient. A polychromatic X-ray source of 100 KeV to 
14 0 KeV is used for penetration of the patient and the 
image is derived from reconstruction within a fast data- 

15 processor of sample data collected from a detector array 
system. The CT contrast number for each voxel of the CT 
image is expressed in terms of the ratio of the 
difference in the X-ray linear attenuation coefficient of 
the tissue from that of the water within the voxel, with 

2 0 respect to that of the water, scaled to 1000. 

In order to optimise the integrated X-ray dose received 
by the patient, and to obtain a reasonable signal-to- 
noise ratio in the image, the sampling is carried out in 
25 accordance with elongate voxels having their longer axis 
normal to the plane of the slice. As a consequence, high 

spatial resolution images are produced only in the 

transverse plane, and reconstruction in oblique or 
otherwise non-transverse planes provides poor spatial 

3 0 resolution . 

A system of radiotherapy-treatment planning (RTP) based 
on CT imaging and called the "EMI Plan", has been used in 
hospitals for many years. However, owing to the 
3 5 inherently poor spatial resolution provided by CT for 

reconstructed images in the non-transverse planes, only 
two-dimensional RTF using CT, namely RTP(CT) , is 



possible. The basic technology used has changed little 
over the years and the limitation on effective imaging to 
allow other than two-dimensional RTP(CT), has remained, 

5 Two-dimensional RTP(CT) has disadvantages which arise 
from limitation of the number of X-ray beam paths that 
can be used, and errors in CT-number computation arising 
from the elongate form of the voxels used for sampling. 
More particularly, there is a theoretical maximum of four 

10 for the number of paths for X-ray beams that can be used 
without creating beam-overlap outside the area to be 
treated; the accumulated dose is required to be as high 
as possible throughout the area to be treated and as 
small as possible in the healthy tissue surrounding, 

15 outside this area. Also, because elongated voxels are 

used for the CT scans, there is a partial-volume effect, 
especially along the edges of a lesion, which in general 
leads to lack of clarity and precision in the 
reconstructed image, impairing visual diagnosis and 

2 0 computation of the isodose curves. 

It is an object of the present invention to provide a 
method of imaging, and an imaging system, of improved 
form as compared with the methods and systems currently 

2 5 used, such as to enable improvement in diagnosis and 

therapy to be achieved. 

According to one of the aspects of the present invention 
there is provided a method of imaging wherein 

3 0 corresponding computed tomography (CT) and magnetic 

resonance (MR) scans of the same part of a subject are 
derived, the scans are related to one another for 
correlation of one to the other positionally with respect 
to said part using a de-convolution mapping process, and 
3 5 imaging of said part of the subject is provided in 

accordance with the MR scan as modified at different 
locations along its length spatially in dependence upon 





the CT contrast numbers applicable to the corresponding, 
correlated locations of the CT scan. 

According to another of the aspects of the present 
5 invention there is provided an imaging system comprising 
means for deriving corresponding computed tomography (CT) 
and magnetic resonance (MR) scans of the same part of a 
subject, means for relating the scans to one. another for 
correlation of one to the other positionally with respect 

10 to said part using a de-convolution mapping process, and 
means for providing imaging of said part of the subject 
in accordance with the MR scan as modified at different 
locations along its length spatially in dependence upon 
the CT contrast numbers applicable to the corresponding, 

15 correlated locations of the CT scan. 

With CT, signal magnitude is monitored by the attenuation 
of the polychromatic X-ray photon energy by all the atoms 
within the beam path. The linear attenuation 
2 0 coefficient, jit, is dependent on the density of the 

absorber and the atomic numbers of the atoms in the X-ray 
beam path, and contrast numbers relative to water on an 
absolute scale are derived. However, MR makes use of 
signals derived from the so-called ^Tl-weighted' and ^T2- 

2 5 weighted' protons of quantum processes, and has contrast 

numbers that are relative in scale. The MR contrast 

numbers are accordingly uncorrelated with those for CT. 

This lack of correlation is confirmed by consideration of 
bone and water which are at opposite ends of the CT 

3 0 contrast (absolute) scale, but which are at the same end 

of the MR-image contrast (relative) scale owing to their 
common low proton-population. 
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The lack of correlation between the CT and MR signals 
acts against their use in combination for imaging 
purposes, but the present invention provides a method and 



system by which the advantages of each may be utilised to 
improve image resolution and contrast information. 

In the latter regard, CT provides a high spatial 
5 resolution but only in regard to view normal to the 
transverse slice. Resolution for all re-constructed 
non-transverse planes is poor owing to the need to use 
elongate voxels to improve signal-to-noise ratio. Also, 
the partial-volume effect of using elongate voxels may 

10 give rise to detection errors at the thin edge of a 

contrast profile of a lesion. MR, on the other hand, can 
give the same high degree of spatial resolution viewed in 
the normal direction to any image-slice plane, and can 
also provide isotropic resolution with cubic-voxel volume 

15 imaging. 




A method of imaging, and an imaging system, in accordance 
with the present invention will now be described, by way 
of example, with reference to the accompanying drawings, 
2 0 in which: 

Figure 1 is a flow diagram illustrating the method, and 
by implication the system also, of the present invention; 



25 Figures 2 and 3 show diagrammatically for the purposes of \ 




explanation of the invention, translation from an object 
domain to an image domain, and a process of convolution 
and de-convolution, respectively; 

30 Figure 4 illustrates relationships between waveforms 

represented respectively at (a) in real space, at (b) in 
image space and at (c) in de-convolution space; 

Figures 5 and 6 illustrate, respectively, one-dimensional 
3 5 and two-dimensional generation of de-convolution maps; 
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Figures 7 to 9 show respectively in perspective, plan and 
side elevation, a form of couch-top used in the method 
and system of the invention; 

5 Figures 10 and 11 show in perspective and cross-section, 
respectively, a form of cylindrical drum phantom for use 
in the method and system of the invention; 

Figures 12 and 13 show in perspective and cross-section, 
10 respectively, a form of elliptical drum phantom for use 
in the method and system of the invention; 

Figures 14 and 15 are cross-sectional views of 
cylindrical and elliptical drum phantoms that may be used 
15 as alternatives to those of Figures 10 to 13; 

Figure 16 is illustrative of signal waveforms derived 
during the method of the invention; and 

2 0 Figure 17 illustrates graphical representation derived 
from the signals of Figure 16. 

In the method of the invention, as indicated in Figure 1, 
multiple-slice transverse CT scans are collected across a 

2 5 section of the volume of interest in a patient or other 

subject. Corresponding multiple-slice transverse MR ^ 

scans of the same volume are collected correspondingly. 

The slice thickness of the latter scans may be one half, 
or smaller, of the thickness of the CT slices, and may be 

3 0 collected two-dimensionally or three-dimensionally . The 

patient is constrained throughout on a couch that 
provides landmarks with respect to a coordinate reference 
arrangement on the couch-top. This is to ensure the 
reproducibility of anatomical positions and features to 
3 5 the first order accuracy for the corresponding CT and MR 
scans, and also for the final radiation treatment to be 
made . 



The respective transverse planes of the CT and MR images 
are processed individually in the method and are matched 
with one another in a de-convoluted space for the CT and 
MR images, as later explained. The two sets of 
5 de-convoluted maps are then merged together to a second 
order of accuracy in order that the CT numbers may be 
transferred over to replace the corresponding MR contrast 
numbers. Once this has been achieved, non-transverse (or 
oblique) planes can be obtained from the two-dimensional 
10 (2D) MR images or from the re-arrangement of the 

corresponding voxels of the three-dimensional (3D) volume 
images; where a 2D-MR is used, a further step of contrast 
transformation may be required. 

15 In CT the transmitted beam intensity I is given by:- 

I = lo exp[-M(x,y)ds] .. (1) 

where: Iq is the incident beam intensity; and 

jLt is the linear attenuation coefficient, which 
is dependent on the density of the absorber 
2 0 (e.g. tissue) and the atomic numbers of the 

atoms in the beam path. 

The integral, ie i:M(x,y)ds, of equation (1) along a ray 
path through the whole body tissue gives the ray sum or 

2 5 the ray projection, p where: 

p = - In (I/Io) • . (2) 

The reconstruction to build up each pixel value in CT 

number is by ^back-projection' and over-lapping of ray 
projections of all the different angles sampled in a data 

3 0 collection. The CT number is given by the expression: 

[{ (jLt (tissue) - M (water)}/ jit (water)] x 1000 (3) 

By adding the effect of the inverse-square law to 
equation (1) above, the percentage dose at depth d below 
3 5 the skin surface, for RTF is given by the expression: 

[{f2 X exp(-Md)}/(f + d)2] x 100 .. (4) 
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where: f is the focus to skin distance (FSD) for the 
photon beam; and 

/i is the linear attenuation coefficient at the 
photon energy of the radiotherapy beam, 

5 

A polychromatic photon energy spectrum is used in both CT 
and RTF for equations (1) and (2) and expression (4) . 

Further in relation to CT, the mass attenuation 
10 coefficient /x/P/ where p is the density of the absorber, 
is given by: 

M/P = Mc/P + T/p + a/p + TT/p (5) 
where: ^i^/p is the mass attenuation coefficient for 
elastic scattering, given by: 

15 Mc/P « zVe . . (6) 

for which Z is the atomic number of the 
absorber and E is the energy of the photon (the 
process involves scattering without absorption 
in the medium, and owing to the low Z in 

20 biological tissues, Mc/P is negligible) ; 

r/p is the mass attenuation coefficient for the 
photoelectric absorption for which: 

r/p oc zVe^ . . (7) 

up to about 2 00 KeV, and at higher energies, 

25 approximates to and eventually to E (in the 

low energy region, the mass attenuation ^ 

coeffigj^nt jg about ^gMiv^lent to th^ mass 

absorption coefficient, and the strong 
dependency on the atomic number of the tissue 

3 0 is the reason for the clear contrast between 

bones and other tissues in X-ray radiography 
and CT imaging - which uses an X-ray beam of 
about 100 KeV - 14 0 KeV) ; 

a/p is the mass attenuation coefficient for 
3 5 Compton scattering and absorption :- 

a/p oc (electron density) /E (8) 
for which electron density = N(A) x Z/A 




N(A) being Avogadro's number, Z atomic number 
and A atomic mass number and: 

ty/p = «y,/p + cr^/p . • (9) 

oJq being the loss by the scattered photon, and 
5 ^J9 absorption loss which is due to the 

energy transmitted to the ^recoiled' electron 
and hence to the medium itself (the energy 
carried away by scattered photons is on average 
greater below 1.5 MeV, whereas the absorption 
10 by ^recoiled' electrons is greater above this) ; 

and 

tt/p is the mass attenuation coefficient for 
pair production: - 

TT/p oc (E-1.02) X Z (10) 
15 for which tt is mainly the absorption of energy 

carried away by two ^annihilation' photons of 
0.511 MeV at the photon energy of 1.02 MeV or 
greater. 

2 0 MR imaging is related to the magnetic moments of nuclei 
(such as the protons of the abundant water molecules and 
tissue proteins present in the body) in precession about 
a uniform magnetic field, BO. This is governed by the 
Larmor Equation: 

2 5 Wo = 7 X BO 

where Wq is the angular frequency of the precessing 
nuclei, and 7 is the magnetogyric ratio of the nuclei. 

The Tl-weighted relaxation relates to the application of 

3 0 radio frequency (RF) to cause the magnetic moments of 

nuclei to be excited to a higher energy quantised state. 
When the RF energy is switched off, the excited nuclei 
return to their lower energy state, releasing the 
absorbed energy; this releasing process, identified as 
3 5 the Tl-weighted relaxation, restores the longitudinal 
magnetisation (spin-lattice interaction between the 
nucleus and its environment) . In practice, the Tl- 




weighted relaxation is measured using an RF field that is 
applied orthogonally to the BO field to rotate the 
precessing nuclei through either 90 degrees (saturation 
recovery) , or more suitably, 180 degrees followed by a 
5 further 90 degrees (inversion recovery) • Detection in 
each case may be from spin echo response or free 
induction decay (FID) • 

The Tl-weighted signal derived from the nucleus in 
10 tissue, is given by: 

(i) using saturation recovery: 

(a) from spin-echo: 

p = p^[l-exp(-{tR-TE}/Tl) ] [exp(-TE/T2) ] (11a) 
15 where: t^ is the repetition time of the successive 

pulse sequences; 

TE is the echo time (from the 90-degree 

pulse to the peak of the spin echo) ; 

Tl is the spin-lattice relaxation time of 

2 0 the nucleus in tissue; 

T2 is the spin-spin relaxation time of the 
nucleus in tissue; 

p is the signal intensity (in arbitrary 
units) ; and 

25 P^ is the value of p at t^ = oo (but in 

practice, t^ may be finite) . ' 

t±g f^T-mn FTD - tmr wh -i r^h is zero reducing 

equation (11a) to: 
p = p^[l-exp(-tR/Tl) ] (lib) 

3 0 or: 

(ii) using inversion recovery 
(a) from spin-echo: 

P = pJl-{2-exp(-{tR-r-TE>/Tl) }exp(-T/Tl) ]exp(-TE/T2) 

. . (12a) 

3 5 where: r is the time between the 180-degree pulse 

and the following 90-degree pulse; 



10 

TE is the echo time (from the 90-degree 
pulse to the peak of the spin echo) equal 
to 2t. 

Tl is the spin-lattice relaxation time of 
5 the nucleus in tissue; 

T2 is the spin-spin relaxation time of the 
nucleus in tissue ; and 

is the value of p at t^ = (but in 
practice, t^ may be finite) • 
10 (b) from FID' - for which TE is zero reducing 

equation (12a) to: 
p = p^Cl-2exp(-T/Tl)+exp{-tR/Tl) ] (12b) 

The second type of signal decay arises from relaxation of 
15 the ^spin-spin' interaction or T2 -weighted signal decay. 

The T2 -weighted proton density and the relaxation in the 
signal decay may be measured by the ^re-focused' or 
^re-clustered' echo-signal or signals. The formation of 
an echo arises if the nuclei are perturbed (or excited) 

2 0 from their equilibrium state (i.e. are rotated away from 

the alignment with the magnetic field, BO) by the 
application of a 90-degree pulse of RF energy. They 
initially all precess in phase with one another (i.e. 
they are coherent in phase) at the switching off of the 

25 RF energy, but with time they are affected by two 

processes, namely: (a) their spin-spin interaction, that 

is, the incremental-field interaction of one nucleus upon 

its neighbours and vice versa, and (b) the inhomogeneity 
of the main magnetic field BO causing the nuclei to 

3 0 precess faster or slower at different locations so that 

phase coherence is lost as time progresses. 

The dephasing effect of (b) can be corrected by applying 
a reverse (18 0-degree) RF pulse after an interval. This 
3 5 will reverse the precessing vectors of the nuclei, thus 
causing the ^re-focusing' or ^re-clustering' of the 
dispersing nuclei to re-establish a coherently-phased 
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signal. The initial transverse magnetisation is thus 
reduced only by the T2 decay, i.e. the spin-spin 
interactions as described in (a) , and can be determined, 



5 The T2-weighted signal magnitude of the nucleus in tissue 
is given by:- 

Pic = P«[l-exp(-{tR-TE(„)}/Tl) ] exp(-TE(k)/T2) (13) 
where: is the signal intensity for the kth of n 

echoes (in arbitrary units) ; 
10 TE is the overall time-of-echo, equal to 2r 

where t is the time between the 9 0-degree and 
180-degree pulses; 

TE(n) is the echo time for the nth echo (from the 
middle of the 9 0-degree pulse to the peak of 
15 the last spin-echo) ; 

TE(k) is the echo time for the kth echo (from the 
middle of the 9 0-degree pulse to the peak of 
the kth echo) ; 

Tl is the spin-lattice relaxation time of the 
20 nucleus in tissue; 

T2 is the spin-spin relaxation time of the 

nucleus in tissue; and 

is proportional to the signal from the 

population of nuclei within a voxel of tissue 
25 being imaged at tR = «> (but in practice, t^ may 

be finite) . ■ 



Other types of scan sequence may be used. 



30 The following TABLE shows the relevant NMR properties of 
various nuclei. 

In the TABLE: RA is percentage relative abundance; 

RS is relative sensitivity at constant 
field for equal number of nuclei; 
3 5 MTR is magnetogyric ratio (7) in MHz per 

Tesla; and 

LF is Larmor frequency in MHz at 1.5T. 
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20 



25 
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TABLE 



Nucleus 








"LF 


1x1 




1 nn 

J. . u 


4 9 58 


63.87 


2H 


0.015 


9.65X10*^ 


6.53 


9.80 


13C 


1.11 


0. 016 


10.71 


16.07 


19F 


100. 00 


0.830 


40.05 


60.08 


23Na 


100. 00 


0. 093 


11.26 


16.89 


3 IP 


100. 00 


6.60x10*2 


17.23 


25.85 


3 9K 


93.1 


5.08x10^ 


1.99 


2.99 




In human tissues, protons from water molecules and 
protein are abundant. The MR imaging technique makes use 
of the signals from the Tl-weighted and T2-weighted 
protons. In the MR scan, the spatial distributions in 
the three dimensions, may be encoded by applying three 
linear gradient fields orthogonal to each other. They 
are generated by coils strategically placed with respect 
to other coils within the system. The RF energy is 
transmitted and received by the transmitter and receiver 
coils strategically placed with respect to the main 
magnetic field, and the gradient fields. The 2D-MR 
sampling uses 90-degree and 180-degree RF pulses for 
excitations with one slice-selection gradient, one 
phase-encode step gradient, one read (frequency encode) 
gradient to obtain a 2D-slice. 



The functions of the gradients may easily be swapped over 
or mixed allowing the data collections for the various 

3 0 orthogonal or oblique slices to be sampled. Hence all 
these slices are equally-well defined in the spatial 
resolution. The three-dimensional volume imaging 
technique uses the non-selective RF hard pulse to excite 
the whole volume, two phase-encode step gradients in two 

3 5 dimensions, with a fixed read (frequency encode) gradient 
in the third dimension. This may provide cubic voxels 
which are isotropic in their spatial resolution. 
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A comparison between CT and MR systems can be summarised 
as follows: 

(a) contrast information in images 
5 (i) CT: Signal magnitude is monitored by the 

attenuation of the polychromatic X-ray photon energy by 
all the atoms within the beam path. The linear 
attenuation coefficient, /x, is dependent on the 
density of the absorber and the atomic numbers of the 
10 atoms. As indicated by expression (3) , the CT number is 
given by the attenuation coefficient for tissue with 
|i reference to that of water to define an absolute scale. 

(ii) MR: Tl- or T2-weighted signal decays of a 
quantum process from a selection of nuclei, e.g. mainly 
15 from protons owing to their relative sensitivity and 

their abundance in tissues (see the Table above) . The 
signal magnitudes are expressed by equations (11) -(13). 
The contrast numbers are relative in scale. 

2 0 By comparing (i) and (ii) , it is quite obvious that the 

contrast numbers for CT and MR are uncorrelated. This is 
clear from consideration of bone and air which are at the 
two opposite extreme ends of the CT contrast scale 
related to water as the base line, but which are at the 

^ 2 5 same end of the MR contrast scale owing to their common 

lack of proton population. Moreover, the CT contrast ' 

number is absolute whereas the MR contrast number is 

relative. 



30 (b) spatial resolution: 

(i) CT: High spatial resolution viewed in the 
normal direction to an image slice plane, is achievable 
only for the transverse slice. All re-constructed, 
non-transverse planes are poor in resolution owing to the 

35 elongated voxel size used in CT because of the 

constraints of minimising accumulated dose and poor 
signal-to-noise ratio. The partial-volume effect of the 
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elongate voxels also gives rise to detection errors at 
the thin edge of a contrast profile of a lesion. 

(ii) MR: The same degree of high spatial resolution 
viewed in the normal direction to any image slice plane 
5 is obtainable by the switching/ swapping of gradient 

fields for 2D-MR imaging. An isotropic resolution of a 
cubic-voxel image may be achieved in 3D-MR imaging. 

(c) motion problems: 

10 (i) CT: Owing to its fast scan, motion blurring in 

the image is arrested. 

(ii) MR: Needs longer scan-time, so is therefore 
more prone to motion-blurring problems in images. When 
motion-suppressed scan sequences or fast scans (such as 

15 the echo planar imaging (EPI) described in GB-A-2235779) 
are developed and used, the motion artefacts can also be 
arrested as in CT; gating may also be used to help. 

(d) geometrical distortion: 

20 (i) CT: This is linear spatially, with no 

geometrical distortion. 

(ii) MR: Geometrical distortions may result from an 

inhomogeneous magnetic field, BO, and/or non-linear 

impurity in one or more of the gradient fields. These 
25 distortions may be corrected by the use of scan sequences 

and/ or phantoms. 

The provision of 3D-RTP has significant advantages. With 
3D there are theoretically up to twelve beam paths for 

3 0 treatment. This greatly improves the absorbed-dose 
ratio, that is to say, the ratio between the dose 
absorbed by the treated area to that absorbed by the 
surrounding healthy tissue, and also provides more 
options for the treatment-beam paths. Accordingly, it is 

3 5 possible to achieve a higher rate of successful treatment 
without serious damage to surrounding tissue or vital 
organs; the possible 300% improvement in the absorbed- 




15 

dose ratio should provide a better treatment success 




The use of MR images in 3D-RTP has the advantage that it 
5 makes high spatial resolution available in all planes, 

and reduces the partial volume effect. Furthermore, the 
use of electronic switching for field variation in MR, 
enables a thinner slice of the 2D orthogonal or oblique 
plane, or cubic-voxel 3D volume, to be achieved. The 
10 combined effect of all the above advantages markedly 
improves the overall diagnosis and treatment success 
rates. 

3D-RTP is not possible using CT in that CT high spatial 
15 resolution is obtainable only in the transverse plane 
and that resolution in reconstructed non-transverse 
planes is poor giving rise to difficulties with both 
diagnosis and the calculation of isodose curves. But, 
and in accordance with the invention, 3D-RTP(MR) becomes 
2 0 a practical possibility because MR provides high spatial 
resolution in all planes. The MR contrast numbers in the 
2D- or 3D-MR images are replaced by the corresponding CT 
contrast numbers, and steps are taken to overcome the 
geometrical-distortion and motion problems. 




The method of the invention, and from this by implication 

the system by which the method is performed, is 

illustrated in the flow diagram of Figure 1, In the 
method, CT multiple-slice transverse (Tr) scans are 
3 0 collected across a section of the volume of interest in a 
patient. Similarly, corresponding MR multiple-slice 
transverse scans of one half, or smaller, of the slice 
thickness of the CT slices are collected across the same 
volume of interest by 2D-MR sampling or 3D-MR sampling. 



The patient is constrained from movement during and 
between the two scanning processes on a universal couch- 
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top which also provides landmarks within a coordinate 
reference arrangement defined in the couch-top. The 
constraint and landmarking is necessary to ensure that 
anatomical positioning and features are maintained the 
5 same to a first order of accuracy for the CT and MR 
scans, and for the final radiation treatment. The 
respective transverse planes of CT and MR images are 
processed individually by using ^ boundary^ or 
^finger-print' matching techniques in a de-convolution 

10 space (DCS) for the CT & MR images. The two sets of 

de-convoluted maps are then merged together to a second 
order of accuracy in order that the CT numbers may be 
transferred over to replace the corresponding MR contrast 
numbers. Once this has been done, non-transverse (or 

15 oblique) planes are obtained from the 2D-MR image or from 
the re-arrangement of the corresponding voxels of the 
3D-MR volume; in the 2D-MR case, a further step of 
contrast transformation is required. 

2 0 Geometrical distortion in the MR images may be corrected 
to a certain extent, by a combined spin-echo (SE) and 
field-echo (FE) scan sequence, with the SE part to 
correct for the field inhomogeneity in BO, and the FE 
part to overcome the impurity or non-linear gradients. 

2 5 The geometrical distortion may also be dealt with 

empirically using a post-process correction based on a ' 
set of phantom data relating to a drum phantom. Movement 
problems may be suppressed by some scan sequences and/or 
by gating; for certain parts of the body there is no 

3 0 movement problem for the MR scan. 




Once the necessary sets of MR images with transferred CT 
numbers have been obtained, 3D-RTP may be started using 
the current 2D-RTP software package. The whole ^3D-' 
3 5 space is covered through the use of several 2D planes in 
a pseudo-3D approach. 
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The results of the combined CT and MR method may be 
accumulated in a data bank so as to enable the technique 
to be performed in a so-called stand-alone MR method. 
For this method the X-tubing coordinate reference 
5 arrangement of the couch top is filled with a chemical 

solution, such as MnClj, of appropriate Tl and T2 values 
to establish a standard in the resulting MR image against 
which the other components of the image can be related. 
In this way, it is possible to standardise the contrast 

10 numbers of various types of tissue using a standardized 
3D-RTP Scan Sequence, A data bank of standardised 
contrast numbers applicable to various tissues of 
different patients versus the CT numbers derived from the 
combined CT and MR method may then be accumulated over a 

15 period of time. The data, as ^banded' in the data base, 
may then be used to enable the stand-alone MR method to 
be performed. Interpretation may still be required by an 
operator or through an ^expert' system to overcome 
certain ambiguities (for example, in separating bone from 

20 air) . 

The matching technique is based on a method using de- 
convolution space (DCS) in combination with a least- 
squares curve-fitting method. These methods are 

25 separately well known in pure mathematics, but not, as in 
the present method and system, combined to produce de- ' 

nnnvolution map s with the least-squares fitting curve 

used as a running filter. The filter is used to overcome 
practical difficulties due to noise within the image 

30 pixels. 

For a mathematical model of the method, consideration is 
given to a line spread function (LSF) arising in the 
imaging system from an impulse or S function; this is 
3 5 represented in Figure 2. 



20 
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An example of the way in which the LSF of Figure 2 is 
convolved with an object step-down function to result in 
a roll-down edge response function (ERF) , and how this 
through the de-convolution process reproduces the LSF, is 
5 illustrated in Figure 3. As illustrated, the LSF is a 

one dimensional profile across the edge of a contour with 
its peak position representing the 50%-level point of the 
ERF. 

10 Figure 3 and the mathematical expressions set out below, 
indicate the way in which the LSF can be recovered from 
the resulting ERF at a contrast edge. For convolution 
related to points 1 to n of the step function: 



15 at point 1: 



at point 2 



J2l{x)Ajc - J^L (x) Ax = E(x)^ 



a-1 



at point 3 : ^ ^ 

5^L(x) Ax - Y^L(x)Ax = E{x)^ 



a-2 



until at point n: ^ 

^Lix) Ax - X)l,(x)Ax = E(x)^ 



25 



For de-convolution ie. for the ERF to manifest the LSF: 

a 

Eix)^ - Eix)^ = 5^ j:(x) Ax = L(x) Ax 



a-1 



thus: 



{E(X)2 - E(X):^)/^X = Lix) (a-2)-(a-i> 



leading to 

iEix) - E(x)J/^x = Lix) 



Thus, the same shape of LSF is recovered by the 
de-convolution process independently of the sense of the 
ERF; it does not matter whether a roll-up or a roll-down 
ERF of a positive or a negative contrast profile, is 
involved. This is an important property as a positive 
contrast contour will have roll-up (from low to high) 
ERFs at both ends whereas a negative contrast contour 
will have roll-down (from high to low) ERFs at both ends. 

In the practical implementation, the LSF is derived from 
the ERF by de-convolution using a running filter. This 
enhances the accuracy of the method in overcoming the 
problem of noise that affects the digitised pixel values 
of the image. Use is made of a least-squares fitting of 
a set of data points that advances along the whole length 
of the function from one point to another. 

Assuming that y = ap + a|X + a,x^ represents the ERF curve 
a five-point or seven-point fit is used, and the normal 
equation becomes: 



El 


Ex, 


Ex? 




aO 


Ey. 




Exf 


Ex] 




al = 


= Ex,y, 


Ex? 


Ex? 


Exl 




a2 


Ex?r, 



25 where: i is 1+n, ... 5+n; ... until m-5, ... m 

or 

i is 1+n, ... 7+n; ... until m-7 , ... m 
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n is 0, 1, 2 , . . • ; and 

1, ... m is the span of the ERF profile, 
The solution may be derived from either: 



aO 


El 


Ex 


Ex2 
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Ey 
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= Ex 


Ex2 


Ex^ 






a2 
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Ex^ 


Ex^ 








or; 



aC7 



Ey Ex Ex^ 
Exy Ex2 Ex^ 
ExV Ex^ Ex^ 



al 



El Ey Ex2 
Ex Exy Ex^ 
Ex2 Ex2y Ex* 



a2 



El Ex Ey 
Ex Ex 2 Exy 
Ex^ Ex^ Ex^y 



.. (14) 



El Ex Ex2 
Ex Ex2 Ex3 
|Ex2 Ex3 Ex 



(15) 



For both equations (14) and (15) to be valid :- 



El Ex Ex2 
Ex Ex2 Ex^ 
Ex2 Ex^ Ex* 



* 0 
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The gradient at dy/dX(3+„) or dy/dX(4+„) can then be derived 
and plotted against x for the LSF profile. 



20 



25 



The graph of dy/dx against x gives the LSF profile. The 
peak o f this pr o file is located half way betwe e n the mid^- 
points of the ascending and descending limbs of the graph 
and is derived from this as the so-called full-width- 
half -maximum (FWHM) with an accuracy to the sub-pixel 
level. 

The point spread function (PSF) is the two dimensional 
profile which may be derived, in practice, from the two 
corresponding LSFs orthogonal to one another within an 
image plane. The peak position of the PSF profile is, 
therefore, from the ^mean' or ^cross-over' of the two 
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peaks or the two LSF profiles. The PSF is obtained in 
practice from two orthogonal axes in a two-dimensional 
plane. 

5 The waveform relationships between real space (RS) , image 
space (IS) , and de-convolution space (DCS) are 
illustrated in Figures 4 at (a) , (b) and (c) 
respectively. These show the relationship between peak 
positions pi, p2 , pi', p2',pl" and p2" and dip positions 

10 d' and d" in the DCS, with the positions and magnitudes 
of contrast value of the waveforms in the image space 
(IS) and real space (RS) , and highlight certain useful 
properties of the process. In particular, the tail-end 
to tail-end distance of each LSF profile gives a spread 

15 of 100% down to 0%, and vice versa, in the corresponding 
ERF profile in IS. Also, the half way points of the 
full-width-half -maximum (FWHM) in the LSF profiles give 
the peak positions of pi, p2 , pi', p2 ' , pi", and p2" 
which correspond to the 50% level points of their related 

2 0 ERF profiles; the dip positions d' and d" in the LSF 

profiles correspond to the central maxima or minima of 
the profiles in IS, and the distances of pi' to p2 ' , and 
pi" to p2" correspond to the FWHMs of the profiles in the 
IS and the distances between the edges of the step 
25 function in RS. 

From the above, the generation of the LSF for PSF) is 

independent of the roll-up or roll-down nature of ERFs at 
the edges of the contrast contour. In other words, it is 

3 0 independent of the sense and the absolute value of the 

contrast numbers within the contour. The peak position 
of the LSF (or PSF) is the central half-way (50%) point 
of the roll-up or roll-down ERF which is the optimum 
position at which the two sets of ERFs are matched 
3 5 together with the overall effect of not distorting the 
contour shape within the boundary of the ERFs. 
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In order to match or merge together the two corresponding 
CT and MR images, they are transformed into DCS, In DCS 
the peaks of the LSFs (or PSFs) mark the central half-way 
points of the ERFs of the original two-dimensional image 
5 planes. These peak positions are then linked together to 
form the ^de-convolution' or ^finger-print' maps. 
Through the assistance of these maps, the two image files 
may now be matched or merged together to a second order 
of accuracy. Figures 5 and 6 illustrate one-dimensional 
10 and two-dimensional generation of de-convolution or 
finger-print maps in DCS. 

The subsequent transformation of the contrast numbers is 
guided by a moving sub-set filter conformed by the 
15 sub-set mean values and the overall uniformity of the 

profile of the MR image; this filter is quite similar to 
the running filter used for the conversion of ERF to LSF. 
The number of samples in the moving sub-set filter may 
also be used for the cut-off rejection of small or 

2 0 ^ point' contours, that is to say, spurious contours 

introduced by noise, within the CT or MR image. 

A similar technique is employed for MR(Tr) to MR (Oblique) 
transformation of contrast numbers. Since a high spatial 
25 resolution and small partial volume effect are 

maintained for the thin 2D-MR (Oblique) slice or the 
3 D-MR (Oblique) plane, the transformation is relatively 
simple. Through this DCS technique, the high spatial 
resolution operation is maintained throughout the whole 

3 0 process of re-processing. 

The combined MR and CT method of the invention requires 
the matching of corresponding CT and MR images with one 
another, and it is important that the same ^laid-out' 
3 5 features of anatomy are presented for the each scanning 
activity, and that the same features are presented again 
for radiation treatment. The couch-top is used for this 
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and the landmarks and reference axes defined in it are 
encoded into the images, for use as guides for the 
repositioning of the patient for any subsequent scanning 
or other activities, and also for the first order 
5 alignment in the initial matching of CT and MR images 
(allowing more refined matching to the second order of 
accuracy using the DCS technique, to be initiated) . 

The form of the couch-top will now be described with 
10 reference to Figures 7 to 9 . 

Referring to Figure 7 to 9, the couch-top has a flat 
upper part 100 that is carried by a curved lower part 101 
which fits to the scanner table of the CT or MR 
15 installation. The patient lies on the flat part 100 and 
is constrained on it (by an arrangement not shown) , so as 
to reproduce the same anatomical features for scanning 
and treatment; the treatment table similarly has a flat 
surface . 

20 

An X-tubing array 102 is embedded into the upper part 
100, The tubing of the array 102 is of plastics or glass 
having an external diameter of up to 6 mm and an internal 
diameter of up to 3 mm, and is filled with a chemical 
25 solution of MnClj* The cross-sectional image of such 

tubing will show a positive profile for MR and a negative 

pr'^fil^ ^'^^ ^"T^ HoWP.Vftr , th<a nr.q tf^nhn-iqng^ nr p-ror!^g.cfi ng 

these images pin-points their centres irrespective of the 
contrast values represented. 

30 

As shown in Figure 8, the patient is constrained by the 
constraint arrangement to lie lengthwise of the flat part 
100 (ie along the z-axis of the coordinate system defined 
in the couch-top) and centrally of it (ie with respect to 
35 the lateral x-axis) ; the patient's head is placed on a 
head-rest 103 suitable for the scan, CT or MR, to be 
undertaken. Three landmarkers/constraint levers 104, two 
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under the armpits and one at the crotch between the 
thighs, are also used to position the patient as well as 
to provide some form of conformation with respect to the 
reference coordinate system. Air-compressed cushions 
5 (not shown) may be gently pressed in from either sides to 

conform the limbs and hence the body along its length. 

The X-tubing array 102 inside the bottom flat part 100 
may be used as the reference system for the transverse 

10 and sagittal planes, and the oblique planes between them. 
The flat surface of the part 100 or the flat contour of 
the patient's skin surface in contact with it may be used 
as the datum for the coronal planes. However, an extra 
small X-tubing array (not shown) may be placed above the 

15 patient using the three landmarkers as the supporting 
legs; this enables the small-angled oblique planes 
between the coronal and transverse planes to be 
monitored. 

20 If the internal-bore diameter of the magnet of the MR 
scanner is large enough, an extra X-tubing array 105 
(Figure 9) may be used along either side of the patient 
to help in monitoring the small-angled oblique planes 
between the coronal and sagittal planes. Furthermore, a 

25 curved X-tubing array (not shown) may be placed inside 
and on the bottom of the head rest (or head coil) 103 

used for the MR scan for the purpose of clearly 

establishing the coordinates of the head. 

3 0 The material of the couch-top is carbon fibre so as to 
ensure a low linear attenuation coefficient with 
transparency to X-rays. Edges of the couch-top are 
smoothly rounded to avoid sharp edges that would cause 
streaking in the CT image. Non-magnetic material is used 

3 5 throughout to avoid disturbing the uniform magnetic 
field, BO, used for MR. 
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The radiation treatment table needs to be capable of both 
translational, "in and out" movement as well as 
rotational movement to cover the 3D planning and 
treatment . 

5 

A final empirical correction of any residual geometrical 
distortion may be carried out using a cylindrical or 
elliptical drum phantom. A cylindrical drum phantom is 
illustrated in Figures 10 and 11, and an elliptical drum 
10 phantom in Figures 12 and 13. 

Referring to Figures 10 and 11, the cylindrical drum 
phantom, which may be of about 160 mm in length and is 
for use with the head coil of the scanner, includes a 

15 central rod 2 00 of polymethyl methacrylate (as sold under 
the Registered Trade Mark PERSPEX) having a diameter of 
about 2 mm and defining the phantom centre, Cph. The 
body of the phantom is built up as concentric sections 
about the rod 200, each comprising evenly-spaced blocks 

20 201 of polymethyl methacrylate and with the blocks 201 of 
each section in exact register radially with the gaps 202 
between the blocks 201 of the preceding section. The 
central region 203, having a diameter of about 20 inm, may 
be left blank (since the magnetic field within the 

2 5 central region is uniform to within 0.1 ppm, so no BO 

correction is required) . The gaps 2 02 and the central ' 

T-ga.g-ion 203 are filled with an NMR solution (e.g. MnCl^ or 

CUSO4) . 

3 0 Peak positions of the PSFs in the DCS map for the drum 

are derived and stored. This data is then corrected for 
each point of the map with respect to the iso-centre of 
the scanner system, Ciso, using the following vector 
equation: - 



35 



Ar = r2 - ri 
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where Ar is the correction vector between Ciso and 

Cph ; and 

r2 and rj are the vectors to the relevant point 
of the map from Ciso and Cph respectively • 

5 

The elliptical drum-phantom of Figures 12 and 13 is 
constructed in essentially the same manner as the 
cylindrical drum phantom, and is for use with the body 
coil of the scanner. The major and minor axes are chosen 
10 according to the field of view being studied. 

As an alternative, the cylindrical and elliptical drum 
phantoms may be constructed as shown in Figures 14 and 
15, using thin plates 3 00 of polymethyl methacrylate 
15 instead of the blocks used for the phantoms of Figures 10 
to 13. The plates 300, which may have a thickness of 
between 1.5 mm and 2.0 mm, are spaced radially from one 
another by gaps 3 01 that, together with the central 
region 302 are filled with appropriate NMR solution. 

20 

CT scan images are not subject to geometrical distortion 
so CT scan data accumulated using the phantoms, and the 
DCS profile data consequently derived and stored, are 
used as a datum or reference set for correction of the MR 

2 5 scans. Such corrections also automatically correct for 
manufacturing tolerances, and the DCS data derived from* 

the phantom through the CT and MR scans represents a 

reference situation devoid of motion. 




3 0 The DCS technique described above is applicable further 
in diagnosis and treatment. In particular, it may be 
used to assist visual detection of gradual change of 
contrast level in scan images. For example, as 
illustrated in Figure 16, the 100% point at which the 

35 roll-off of a lesion profile begins can be readily 

plotted from the DCS profile, together with the 50%-level 
point, and then the 0% point where the profile has 
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completely levelled off. Graphic representations of 
these various points can then be plotted out to 
facilitate diagnosis, as illustrated in Figure 17, in the 
form of a string of characters 400 delineating the 100%- 
5 level, a string of characters 401 delineating the 50%- 
level and a string of characters 4 02 delineating the 0% 
level . 

In current radiation treatment planning, a composite 
10 volume leaving a margin of 15 mm around a lesion, is 

used. To be effective, the tail-end of the lesion growth 
is to be within the treatment-beam path, and in this 
regard there is advantage in extending the DCS-profile 
plot of Figure 16 to include a point which is a distance 
15 (d) beyond the 0%-level. This is then used for diagnosis 
and treatment, to establish an outer string of characters 
403 in the graphical representation of Figure 17. The 
distance (d) covers the penumbra-spread as well as any 
possible mis-alignment by the operator. 
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